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ABSTRACT 


Effect of Soil-structure interaction (SSI) modelling on seismic vulnerability of reinforced concrete moment resisting frame building is 
the key objective of the present study. The building situated in Indian seismic zone IV and resting on soft soil has been considered 
for the study. Three types of foundations have been assigned to example building; fixed base, isolated footing and raft foundation. 
For SSI effect for the building with isolated footing and raft foundation, a linear spring modelling have been used. Seismic inelastic 
displacements in these three models have been estimated from static pushover analysis. Spectral displacement based fragility curves 
have been developed for different damage states. Seismic reaction of building considering SSI modelling is related with 
conventional fixed base model. It is concluded that for the building resting on soft soil, SSI modelling is the influencing parameter 
on base shear, drift, ductility and on probability of exceeding the damage state. 
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1. INTRODUCTION 


Seismic fragility curves are the tool to asses seismic vulnerability of existing buildings and several researchers were used it for the 


estimation of losses due to earthquake (Nazri, 2018). These curves are generated from real earthquake damage data to estimate or 
predict whether the damage meets or exceeds a certain performance level for a given set of ground motion parameters (Pitikalis et 
al. 2014). In addition, the curves can be applied to predict both pre- and post-earthquake situations. These curves are unique 
because every building has specific fragility analysis (Hancilar et al. 2014). 

Interactions between three linked systems: the structure, the foundation, and the soil underlying and surrounding the foundation 
are the reason to affect the response of the structure while ground waves excitation (FEMA440, 2009).Soil-structure interaction (SSI) 
study estimates comprises the effect of a specified free-field ground motion, on the structure, the foundation and the soil conditions 
around the foundation. 

SSI effects are absent for the theoretical condition of a rigid foundation supported on rigid soil. Accordingly, SSI accounts for the 
difference between the actual response of the structure and the response of the theoretical, rigid base condition (NIST GCR 12-917- 
21, 2012). 

This paper focuses on SSI analysis considering various foundations. The influence of SSI modelling on seismic vulnerability of 
structure is studied by evaluating fragility curve. 


2. ANALYTICAL METHODOLOGY 


Menglin et al. (2011) concluded from the research that SSI is beneficial to explore seismic response in true sense. The main focus 
was to evaluate the effects on SSI on the structural response with consideration of seismic provisions. Mylonakis & Gazetas (2000) 
concluded that SSI affects on ductility demand of the structure which may either increase or decrease depending on the features of 
ground motion and structure. 

For present study, analytical method based on static pushover is adopted to develop fragility curves. A static pushover method has 
been used to estimate seismic inelastic displacements in the structure. The direct analysis approach is used for modelling soil in 
which soil and structure are included in the same model. Winkler model representation is used to assess the effect of soil lying 
underneath the foundations. Linear springs with six DOF (three translations and three rotations) were assigned for modelling the soil 
(Saez- 2012, ASCE 41-2017). A parametric study to assess the seismic response of building has been carried out by considering 
different types of foundations resting on soft soil. 


36 36m 
m 


(height) 


17m 
(Length in X) 17m 


Figure 1. Elevation of model a) building with fixed base, b) building with spring supports 
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Example Building Model 


A 3-D frame building having five stories has been considered. The elevation of it is shown in Figure 1. It is an industrial building 
having storey height as 5.5m each. An importance factor of 1.5 as per IS 1893-2016 and the building measuring 17 m X 40 m in plan 
has been used. Figure 2 show the plan of the building. The building is designed as a special moment resisting frame using linear 
response spectrum analysis and designed for various load combinations specified in the code. Concrete of grade M30 and 
reinforcement having grade Fe500 have been considered. Table 2 shows the details of column sections. The beam sections are 
230X450, 230X600, 230X750, 300X400, 300X500, 300X600, 300X700, 300X750 and 450X750 mm. 

Three models have been considered for the study a) Model I- model with fixed base, without SSI b) Model Il- model considering 
SSI with isolated footing at base c) Model IIl- model considering SSI with raft foundation at base. 


40 m 


17m 


Figure 2. Plan of the building model 


Modelling of foundation 
Sizes of the isolated footings are mentioned in Table 2. The embedment depth has been assumed to be 2m.Thickness of raft 
foundation is 200mm. The equations from ASCE 41-2017 have been used to calculated equivalent spring stiffness’s. Soil properties 
in terms of the elastic modulus, shear modulus and Poisson’s ratio were considered appropriately in the calculation of equivalent 
spring stiffness. Table 1gives the soil properties considered for modelling. Table 3 shows the soil spring stiffness values used for 
both Model II and Ill. 

For modelling of isolated footings in Model Il, six springs were assigned at the end of the columns as shown in fig. 1b. Raft slab 
in Model Ill is modelled as area- spring modelling command in SAP2000v21 software (Raychowdhury, 2008). 


Table 1. Properties of soil 


Type of soil Soft soil 

Soil bearing capacity 100 MPa 
Young's modulus 12000 MPa 
Poisson’s ratio 0.45 

Shear modulus 4137.93 MPa 


Table 2. Reinforcement details in Columns 


Column Size in mm Main reinforcement Footing size in mm 
C1 500X500 12#20o 2250x2250x750 
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C2 550X550 12#25 2300x2300x750 
C3 600X400 6#25 b and 6# 20 > 2600x2600x850 
C4 650X450 12#25 2700x2700x850 


Table 3. Spring stiffness values at foundation 


Spring stiffness at foundation Values(KN/mm/) 
Translation along X direction - Kx 160745.4 
Translation along Y direction - Ky 165040.4 
Translation along Z direction - Kz 2329174 
Rotation along X direction - Rx 2668981 
Rotation along Y direction - Ry 4202134 
Rotation along Z direction - Rz 3721987 


3. MODAL ANALYSIS 

Modal analysis has been performed to study the dynamic properties of the building, Table 4 shows the result. It has been observed 
that the time period for model | is less than other models, indicating stiff behaviour resulting conservative design. Mass participation 
in fundamental mode for all models is less indicating need of higher modes for response estimation. The coupling of modes has 


been observed. Analysis results are also showing that X direction is weaker. 


Table 4. Modal Analysis Results 


. ; Mass participation in fundamental 
Fundamental time period : 
Model (Sec) mode (unit less) 
UX Rx 
Model | 2.04 0.73892 0.10523 
Model Il 2.14 0.74803 0.10472 
Model Ill 2.13 0.6347 0.1495 


Nonlinear modeling of structural elements 

As per to structural mechanism concept, ends of beam and column are the likely locations of formation of hinges under lateral loads 
due to concentration of flexural and shear stresses. Flexural and shear failure of the elements are commonly observed; they are 
usually addressed by assigning plastic hinges, near the ends of the element, one is flexural hinge and other is shear hinge. 

As the building considered under study is modelled according to 1S1893 and IS 13920 revised code; it is safe guarded for shear 
and hence there is no need to assign shear hinge for the beam and column sections as well as joint section. Flexural hinges of M3 
and P-M2-M3 DOF have been assigned to beams and columns respectively. In SAP2000, these hinges can be assigned directly using 
auto-hinge option. Slab has been assigned with rigid diaphragm. Geometric nonlinearity has been considered by incorporating P- 
delta effect. Contribution of brick masonry infill in resisting lateral loads has neglected. Nonlinear Pushover analysis (POA) has been 
performed to develop seismic fragility curves. A capacity spectrum method has been used to estimate seismic performance and 


vulnerability of the structure. 


4. PUSHOVER CURVES 

Pushover curves for the model |, Il and model Ill were obtained from POA. These curves along X and Y directions are shown in 
Figures 3 and 4. The base shear of the structure is affected by different foundations considered. Base shear for the model | is 
comparatively higher than other two models indicating higher lateral resisting capacity. The roof displacement of the building 
increases by 11% for model Il and III. Thus by consideration of SSI, the behaviour alters from stiff to flexible. 


Table 5. Base shear at performance point 


Model Base Shear (KN) 
X- Direction Y-Direction 
Model | 7511 5770.162 
Model II 6905 5524.618 
Model III 7143 5593.397 
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PUSH OVER X CURVE 
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Figure 3. Push over curve in X- direction 
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Figure 4. Push over curve in Y- direction 


Drift and Ductility 
Based on the POA results, inter storey drift are calculated for model |, model Il and model Ill. Corresponding graphs in X and in Y 
directions are as shown in Figures 5 and 6 respectively. 

It has been seen that the storey drift increases for model Il and model Ill by 6 times. Model Ill having raft footing has 
comparatively less drift than model II. It shows that in soft soil conditions, isolated footing (model II) is more susceptible to damage 
then raft footing (model III). Pushover curve is converted into bilinear curve (ATC40-1996). The ductility factor of model |, Il and Ill 
were calculated from bilinear curve as a ratio of ultimate displacement to the yield displacement. Table 6 shows the result. The 
ductility of structure increases in model II and model Ill due to consideration of SSI effects. 
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Table 6.Ductility factor 


Model Ductility factor in X-Direction (Unit 
less) 
Model | 10.96 
Model II 12.20 
Model Ill 11.83 
STOREY DRIFTIN X 
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Figure 5.Storey drifts in X- direction 
STOREY DRIFTINY 
—MODELI —MODELII ——MODELIII 
8 
7 
6 
Be, 
O4 
ta 
an 3 
2 
1 
0 
0 0.0005 0.001 0.0015 0.002 0.0025 0.003 0.0035 
STOREY DRIFT 
(UNITLESS) 


Figure 6.Storey drifts in Y- direction 


Damage Pattern 

The damage pattern has been studied for three models by POA at performance point, number of hinges developed in the structure 
is shown in Table 7. At performance point, the number of hinges developed in model | are in advanced state than model II and Ill. 
Thus the assumption of fixed base shows conservative performance. 
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Table 7. Hinge states at performance point 


Model Direction Number of hinges 

ATOB BTO IO 10 TO LS LS TO CP CPTOC Total 

Model | PUSH X 596 550 778 590 0 2514 
PUSH Y 567 232 934 781 0 2514 

Model II PUSH X 1779 409 217 109 0) 2514 
PUSH Y 1784 207 310 213 0) 2514 

Model Ill PUSH X 1409 458 416 231 0 2514 
PUSH Y 1397 201 577 339 0 2514 


5. FRAGILITY ANALYSIS 


Fragility curves are defined as the probability of the exceedance of the given damage states d,j, for the spectral displacement of the 
structure, which ultimately shows the vulnerability of the structure. 


A fragility curve is well described by lognormal probability density function shown by eq.1 (HAZUS MH MR4 -2003). 


. (1) 
P(dsi/Sa) = 9 hie 7 (=| 


Where Sd,iis the threshold spectral displacement at which the probability of the damage state ds; is 50 %, Basi is the standard 
deviation of the natural logarithm of this spectral displacement, @ is the standard normal cumulative distribution function. The 
damage states threshold is calculated using the equations shown in Table 8. 

The fragility curves were plotted using the data from capacity spectrum analysis. The pushover curve is idealized into bilinear 
curve as per FEMA 356-2000 coefficient method. The spectral displacement at yield and ultimate point were obtained from bi-linear 
curve. 


Damage state definition 

Building damage functions are in the form of lognormal fragility curves that relate the probability of being in or exceeding. Four 
damage states (DS) have been considered in present study; slight damage (SD), moderate damage (MD), extensive damage (ED) and 
complete damage (CD). Damage state thresholds (Sq,;) for each DS, median threshold value of Sa has been estimated from Figure 


7. Table 8 shows estimated median values for each state. 


Spectral acceleration, Sa 


Sd, Sd, 3d, Sd, 
Spectral displacement, Sd 


Fig 7 Damage State thresholds for capacity spectrum (Barbat et al 2008) 


Table 8 Damage state thresholds (S,,;) 


Damage state thresholds 
Sdi= 0.7Dy 
Sd> = Dy 


Damage state 


Slight 
Moderate 
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Sd3 = Dy+0.25(Du-Dy) Severe 
Sd4 = Du Complete 


Estimation of fragility curve parameters 
Fragility analysis has been carried out for weaker direction (X-direction). Median spectral displacement (Sq,;) estimated values are 
shown in Table 9.The total variability (Bg.;) value is taken from HAZUS (2003) manual. 


Table 9. Estimated parameters of fragility curves 


Median parameter(S,,,) (cm) Variability parameter 
Damage states 
Model | Model II Model III (Basi)(cm) 
X-direction X-direction X-direction For X direction 
Slight 0.0518 0.055386 0.05458 1.6764 
Moderate 0.074 0.079123 0.077972 1.6256 
Extensive 0.2584 0.29317 0.295157 1.7018 
Complete 0.8117 0.93523 0.93972 1.9812 


Fragility curves and discussion 
The fragility curves were developed for three models. Figures 8, 9 and 10 shows the fragility curves for model I, Il and Ill respectively. 
The Table 10 shows the comparison of probability of exceedance for collapse damage state at 2 cm spectral displacement. 


Table 10. Probability of exceedance at 2 cm Spectral displacement 


Probability of exceedance (%) 
Model SD MD ED cD 
X-Direction X-Direction X-Direction X-Direction 
Model | 98.53% 97.8 % 88 % 67 % 
Model II 98.38 % 97.6 % 87 % 64.9 % 
Model Ill 98.4 % 97.7 % 86.95 % 65 % 


FRAGILITY CURVE IN X FOR MODELI 
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Figure 8. Fragility Curve in X - Direction for Model | 
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FRAGILITY CURVE IN X FOR MODELII 
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Figure 9. Fragility Curve in X - Direction for Model Il 
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Figure 10. Fragility Curve in X - Direction for Model III 


6. CONCLUSION 


Due to SSI following effects are seen: 


a) 
b) 
Cc) 
d) 
e) 


Roof displacement increases by 6 times for models with SSI effect than fixed base model. 
Base shear decreases by 4.8 times for models with SSI effect than fixed base model. 

Time period increases by 4.9 times for models with SSI effect than fixed base model. 
Inter-storey drift increases by 6 times for models with SSI effect than fixed base model. 
Ductility increases by 7.9 times for models with SSI effect than fixed base model. 
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f) Probability of exceedance for fixed based model is higher in ED and CD state than models with SSI effect 


Concluding remark 

This study proved the importance of consideration of soil-structure interaction on performance-based design of buildings and the 
raft foundation gives quite similar results to fixed base rather than isolated footing. Economic design can be achieved by including 
SSI effects in seismic design. Study recommends SSI modelling for the structure resting on soft soil since it is realistic case. 
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